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Abstract 

We calculate the light-cone wavefunctions of the pion by solving the meson boundstate problem 
in a coarse transverse lattice gauge theory using DLCQ. A large- A'^c approximation is made and the 
light-cone Hamiltonian expanded in massive dynamical fields at fixed lattice spacing. In contrast 
to earlier calculations, we include contributions from states containing many gluonic link-fields 
between the quarks. The Hamiltonian is renormalised by a combination of covariance conditions 
on boundstates and fitting the physical masses Mp and Aij^, decay constant /^r, and the string 
tension y/a. Good covariance is obtained for the lightest state, which we identify with the 
pion. Many observables can be deduced from its light-cone wavefunctions. After perturbative 
evolution, the quark valence structure function is found to be consistent with the experimental 
structure function deduced from Drell-Yan pi-nucleon data in the valence region x > 0.5. In 
addition, the pion distribution amplitude is consistent with the experimental distribution deduced 
from the Trj*j transition form factor and diffractive dissociation. A new observable we calculate 
is the probability for quark helicity correlation. We find a 45% probability that the valence-quark 
helicities are aligned in the pion. 
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I. INTRODUCTION 



Light-cone wavefunctions encode all of the bound state quark and gluonic properties of 
hadrons, including their momentum, spin and flavour correlations, in the form of universal 
process- and frame-independent amplitudes (see, for example, Ref. [jll). Hadronic observ- 
ables represented as matrix elements of currents are easily expressed in terms of overlaps of 
light-cone wavefunctions. To compute the wavefunctions, one must diagonalise the light-cone 
Hamiltonian of QCD in a Fock space of quark and gluonic degrees of freedom. A promising 
method to achieve this is the transverse lattice formulation of gauge theory 0, ^. In this 
approach, the physical gluonic degrees of freedom are represented by gauge-covariant links 
of colour flux on a lattice transverse to the null-plane of quantisation. In this paper, we set 
up the method and solve for the light-cone wavefunctions of light mesons using a physically 
realistic truncation of Fock space on a coarse lattice, spacing ~ 2/3 fm. We obtain good 
covariance for the light-cone wavefunction of the lightest meson, which we identify with 
the pion. Results for the pion distribution amplitude (valence quark wavefunction at small 
transverse separation) and distribution function (valence quark probability at any transverse 
separation) are consistent with the most recent experimental results in the valence region of 
light-cone momenta. We find the distribution amplitude to be 

(f)^{x) = 6x(l -x)[l + 0.15(2) 2(1 - 2x) + 0.04(1) Cf ^(1 - 2x)} , (1) 

while the distribution function is 

^-(^) = ^0.7(1) {0-33(3) - 1.1(2)7^ + 2.0(3)a;} . (2) 

where x is the quark light-cone momentum fraction carried in the pion. The transverse 
renormalisation scale should be 0.5 GeV if the first moment of K- is to agree with experiment. 
As a further application of the light-cone wavefunctions, we also compute the probability 
for a valence quark of momentum fraction x to have its helicity correlated with that of the 
anti-quark in the pion. We find a surprisingly large probability ~ 45 % for the quark and 
anti-quark helicities to be aligned, even though the pion is spin 0. These represent our main 
results. 

Attempts to solve transverse lattice QCD have been renewed in recent years for both 
the pure gauge theory |^, ^ and mesons 0, |], |^. The most succesful approaches have 
employed the original idea of a l/A^^c and colour-dielectric expansion in dynamical fields 
to approximate the light-cone QCD Hamiltonian on a coarse transverse lattice. For pure 
gauge theory, to lowest non-trivial order of the expansion, requirements of vacuum stability, 
Lorentz and gauge invariance alone were found to constrain the coarse lattice Hamiltonian 
sufficiently accurately for first-principles predictions of the glueball states Extension of 
this work to light mesons introduced quarks and imposed a (Tamm-Dancoff) restriction on 
the number of link fields in Fock space In previous calculations ^], no more than one 
link field was allowed in a meson. This effectively restricts the transverse size to < 2/3 fm, 
which is unrealistic for light mesons. In this case, the correct Hamtiltonian could not be 
accurately identified using Lorentz and gauge invariance alone. Some phenomenology was 
also needed. 

In this paper, we again use the lowest non-trivial order of the colour-dielectric expansion 
of the hamiltonian, but relax the Tamm-Dancoff cut-off on the space of states. This allows 
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light mesons to expand to their physical transverse size. It also means that one begins to 
take account of the full pure-glue dynamics in the meson sector. While the results are now 
realistic, we find that it is still necessary to use some phenomenological fitting of masses 
and decay constants, in addition to optimizing Lorentz covariance, to obtain unambiguous 
couplings in the coarse-lattice Hamiltonian. We believe this is due to the absence, in the 
currently employed transverse lattice Hamiltonian, of operators needed to optimize chiral 
symmetry. We show that such operators would occur at higher order of the colour-dielectric 
expansion. In the next section, we review and extend the previous work. Section |ITT| describes 
the procedure we employ for fixing the various couplings that appear in the Hamiltonian. 
Finally, our results for pion observables are discussed in Section |V|. Chiral symmetry issues 
are discussed in an Appendix. 

II. TRANSVERSE LATTICE MESONS 
A. Hamiltonian 

We introduce continuum light-cone co-ordinates = (x" ± x^)/-\/2 and discretize the 
transverse coordinates x = on a square lattice of spacing a. Lorentz indices /z, v are 

split into light-cone indices a,/? G {+, — } and transverse indices r, s G {1,2}. Subsequent 
analysis is done to leading order of the 1/Nc expansion of the gauge group SU{Nc). Quark 
fields \l/(x+,x~,x) in the fundamental representation and gauge potentials y4a(x+,x~,x) in 
the adjoint representation of SU{Nc) are associated to the sites of the transverse lattice. 
Link fields Mr{x^ ,x~ which we choose to be complex N^xNc matrices, are associated 
with the directed link from x to x -|- of. They carry flux from site to site. This use of 
disordered link variables implies that a coarse transverse lattice is being considered. 

For finite spacing a, the Lagrangian can contain any operators that are local, invariant 
under transverse lattice gauge symmetries and under Poincare symmetries manifestly pre- 
served by the lattice cut-off, and renormalisable by dimensional counting with respect to 
the continuum co-ordinates x°'. The objective is to obtain an approximation to the light- 
cone Hamiltonian operator P~ , that may be diagonalised in a Fock state basis of the above 
fields. This may be achieved by first fixing to the light-cone gauge A_ = 0, eliminating 
non-dynamical fields, then expanding the resulting Hamiltonian in powers of the remaining 
dynamical fields. Truncation of such a 'colour-dielectric' expansion is a valid approximation 
provided wavefunctions of interest (typically those of the lightest eigenstates) are dominated 
by few-body Fock states. This is achieved by working in a region of coupling space with 
sufficiently heavy dynamical fields. This in turn will be found to constrain the transverse 
lattice spacing a to be coarse when symmetries and phenomenology are optimized. 

The Lagrangian density we consider contains terms up to order and \1/M\1/ for the 
large- iVc theory, 

L = T.fdx- E E -^^TrfF-'^lxjF^.lx)} 

X •' a,l3=+,- r=l,2 

+ Tr{D,M,(x)(D"M,,(x))t} 
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(^'(x)7''M^(x)^'(x + ar) - ^(x)7''M/(x - ar)^(x - ar) 
+Ks (*'(x)M^(x)^(x + ar) + "^(x)M^t(x - ar)^(x - ar)) - , (3) 

where F°^(x) is the continuum field strength in the {x^, x^) planes at each x, 

D„M,(x) = {da + iA«(x))M,(x) - iM,(x)A^(x + ar) (4) 
and the link-field potential is 

K = --p^ETr{M,(x)M,(x + ar)Mt(x + as)Mt(x)} 

r^s 

' Tr{M,MtM,Mt} + ^ Tr{M,(x)M,(x + ar)Mt(x + ar)Mt(x)} 



+7^ E Tr{MtM.Mi,M.4 . (5) 

We have defined Mr = mI^, and hold G G^/Wl finite as Nc — >■ oo. To this action we 
could in principle add allowed operators at order M^, (\l/\l/)2, ^^M^vl/, and so on. It should 
therefore be understood as the truncation of an expansion in powers of the fields. Strictly 
speaking, this expansion should be performed for the light-cone gauge-fixed Hamiltonian in 
terms of dynamical fields only. 

In the chiral representation, = ("U^j., fi^, f* , m1)/2^/^ decomposes into complex fermion 
fields V {u) with a helicity subscript /i = ± denoting the sign of the eigenvalue of 7^. In 
light-cone gauge A_ = 0, and v± are non-dynamical (independent of light-cone time x~^) 
and are eliminated at the classical level using the equations of motion 

id-)'A^ = ^{-J^-j^T^'J^) (6) 
id-v, = ^F., , (7) 

where we have defined 

F/j(x) = -M/,(x) + — ^ (^Mr(x)M/i(x + ar) + A4(x - ar)M/,(x - ar)j 

(Mi(x)M_ft,(x + ai) - hiM2{x))u_hi^ + a2) 
-MI{x - ai)u_hix - ai) + /iiM|(x - a2)M_ft(x - a2)} , (8) 
J+(x) = i^ ('m,(x) d- Mt(x) + Mt(x - ar) d- M,(x - ar) 

+ ^w,(x)4(x), (9) 

h 

The lightcone Hamiltonian, expressed in terms of the remaining dynamical fields u±(x) and 
Mr(x), may be obtained from the action (^) in the standard way 

+yE(^/I^^^)+^x[m]. (10) 
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Under certain reasonable assumptions [§, the Hamiltonian ([T0| ) is a truncation of the 
most general Hamiltonian to order M*^ and uMu. It also contains some, but not all, al- 
lowed operators at order uM'^u and m^. In particular, it contains the combination J^dz^J^, 
which is responsible for confinement in the lattice theory of states singlet under residual x~- 
independent gauge transformation 0. The various parameters G, ^f, Kq, k^, fit,, Ai, A2, A4, j3, 
as well as ones that would appear at higher orders of the colour-dielectric expansion, are 
coupling constants that will run with the cut-off(s) in the theory. In principle, this running 
could be determined by performing renormalisation group transformations from QCD at 
short distance scales. However, on a coarse lattice, weak-coupling perturbation theory is 
not available, and such an approach become unworkable. One may also treat the problem 
as an effective field theory, fixing couplings phenomenologically. Even in this case, one may 
constrain the parameters from first principles by empirically tuning them to minimize the 
violation of continuum symmetries. In the case of pure gauge theories, at lowest order of 
the colour-dielectric expansion, this gave a quite accurate estimate of the running couplings, 
without the need to resort to 'phenomenology' For meson calculations with our choice 
of Hamiltonian (p!OD, additional phenomenological constraints must be used to obtain un- 
ambiguous values for the coupling constants, although symmetry requirements do strongly 
constrain them. 

Of the other generators of the Poincare algebra, P'^,M'^'', the following can be derived 
canonically at x"*" = 

P+ = y" dx~ J2 2Tr{9_M,(x)9_M,(x)"^} +i<9_M/, , (11) 

x,s,h 

M-+ = j dx-Y. X' |2Tr{9_M,(x)9_M,(x)t} + ^< , (12) 

-x.,s,h 

= - Jdx-Y'^(^' + -'^") Tr{9_M,(x)a_M,(x)t} + ix'uld.Uh . (13) 

x,s,h 

Note that these are all kinematic operators, quadratic in fields. and M generate 
translations and boosts respectively in the x~ direction and are unaffected by the transverse 
lattice cut-off. The cut-off effects on the boost-rotation operator M~^^ are discussed further 
in the next section. 



B. Space of states 



For the construction of a Fock space of the dynamical fields Mr and Uh, it is convenient 
to Fourier transform the fields in the x~ co-ordinate only. Thus, we introduce a Fock space 
operator al^j{k~^ ,x.) which creates a 'link-parton' with light-cone momentum k~^, carrying 

colour i G {1, ■ ■ ■ , Nc} at site x to colour j at site x + ar; aL 

K 

at site x, while d* does the same for anti-quarks. We have 



_^ creates an oppositely oriented 
link-parton. Likewise, b'^^{k~^,x.) creates a 'quark-parton' of helicity h, colour i, momentum 



ax,ij{k~^,:>i.),a*^i{k^,y = 6ikSji6xpS^y6{k'^ - k~^ 



ax,ijik^,x),apMik~^ 



0, 



(14) 
(15) 
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{bh,{k^,^),bl,/k+,y)} = 5i,5Hh'S^y5ik+-~k+) , (16) 
{6;,,(A:+,x),6^,,(^+,y)} = 0, (17) 

where A and p G {±1, ±2} denote the orientations of hnk variables in the (a;^,a;^) plane, 
~ ^"^^ similar anti-commutators exist for d. Fock space is already diagonal in 

the light-cone momentum P"*" and serves as a basis for finding the eigenfunctions P~ , the 
light-cone wavefunctions. As usual in light-cone quantisation (without zero modes), the 
Fock vacuum state |0) is an exact eigenstate of . 

Further cut-offs, apart from the transverse lattice, must be applied to Fock space to make 
it finite-dimensional. We will use DLCQ |Tl| to discretize light-cone momentum, which 
amounts to compactifying x~ on circle of circumference C = 27iK/P~^, where K is a positive 
integer, with periodic (anti-periodic) boundary conditions for M {u). Eventually, we will 
extrapolate observables to = oo. The use of anti-periodic boundary conditions is desirable 
because it tends to improve convergence a.s K oo. However, one cannot consistently have 
anti-periodic boundary conditions for both bosons and fermions in a theory with Yukawa- 
type interactions. 

To reduce the size of Fock space still further, it will be convenient to impose a separate 
Tamm-Dancoff cut-off on the maximum number of partons in Fock space, studying the theory 
as this cutoff is raised. Since the large limit automatically restricts to a quark- ant iquark 
pair in the meson sector, this effectively means a cut-off on the number of link-partons. A 
general meson state of light-cone momentum P+, which is translationally invariant in the 
transverse direction, takes the form 

|^(P+)) = 

2av/?^^ /^''^ciA;+d4^(P+-A;+-A;+){^w(a:i,a:2) 6l(A;+, x)4,(fc+, x)|0)} 



/N 

^^^c X h,h' 

E E ^^^|^5(p+ - kt - kt - kt 



)/j/(xi,X2,X3) 6[(A;]*",x)a|;(/i;^,x)rf^,(A;^,x + ar)|0) 
+ iJh{-r)h'{xi,X2,x^) 6|^(A;^,x + ar)a;Lr(fc^,x)(i^,(fc-[,x)|0)} 
2a^ ^ ^ fP^ dktdktdk^dkt 



+^ EE/ " 7;: " " siP^ - K - kt - kt - kt) (18) 



rs)h'{xi, X2, xs, X4)bl{kt, x)aJ(A;^, x)aJ(A;^, x + ar)dl,{kt, x + ar + as)|0) 
+iph{r-s)h'{xi, X2, xs, X4) b[{kt,yi)al{kt,^)a'Ls{kt,y:^ + ar - as)dl,{kt ,y:^ + ar - as)\0) 
+A{-rs)h'{xi,X2, X3, 0:4) blikf, X + ar)al^{kt, x)al{kt, x)dl,{kt, x + as) |0) 
+il)h{-r-s)h'{xi,X2,X3,X4,) 6j,(/c+,x + ar)air{kt , ^)a-s{kt , ^ - as)(i^,(A;|, x - as)|0)| + 



where states with up to two links have been explicitly displayed. In (^Sf ), f acts on gauge 
indices and Xi = kf/P^ etc., are light-cone momentum fractions. Only gauge singlet 
combinations under residual gauge transformations in A_ = gauge can contribute to 
states of finite energy 0. Because pair production of quarks and mixing with glueballs is 
suppressed at large N^, the states (|^) provide a description of the valence quark content 
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of flavour non-singlet mesons. Thus one should implicitly understand a distinct flavour 
label on the quark and anti-quark, which is redundant. For states that are translationally 
invariant in the transverse direction, the transverse co-ordinate argument in wavef unctions 
may be suppressed. The sequence of orientations A, p, ■ ■ ■ of link variables, together with 
the P"*" momentum fractions xi,X2, ■ ■ ■ and quark helicities h, h' are sufficient to encode the 
structure of Fock states contributing to the boundstate. Thus, a general Fock state may be 
labelled 

|(xi,/i),(x2,A),---,(x„_i,p),(x„,/i')) • (19) 

The expansion (|1^) may be represented by a planar (large- A^c) diagram notation shown in 
Fig. |l]. This will be helpful when enumerating the matrix elements of the Hamiltonian. 

The transverse momentum operator is not directly defined because of the lattice regulator, 
but one may introduce transverse momentum P by application of the boost-rotation operator 
M+^. Let |(xi, xi), ■ ■ ■ , (x„, x„)) denote an n-parton Fock state. Xp is the transverse position 
and Xp the momentum fraction of the p^^ parton (conventionally we take transverse 
position to be the midpoint of a link, for link fields). Using (|13D we find 



exp 



iM+^'P,/P+l |(xi, xi), ■ ■ ■ , (x„, x„)) = exp 



iP. ^ ^ XpX.p 



(xi,Xi), ■ • - , (x„,x„)) . 



(20) 

Therefore, the net effect is to add phase factors into matrix elements of P between Fock 
states at P = 0. In a Poincare covariant or a free theory, the transformation (^) applied 
to eigenstates of P^ (|l^) would be sufficient to generate eigenstates of P~ at non-zero P. 
However, the lattice regulator spoils Poincare covariance and in general one must rediago- 
nalise P~ after boosting Fock states by (pO|). Thus the eigenfunctions if) in (pIS]) for P~ will 
become functions of P also. 

The state is normalised covariantly 

(V^(P+, Pi)|V^(P2+, P2)) = 2P+(27r)35(P+ - P2+)5(Pi - P2) (21) 

if 

1=1 dxY\i^hh'{x,l - X)\^ + I dXidX2 Y \i^h(X)h'{Xi,X2,l - Xl - X2)\^ 
•'^ h,h' ''^ h,X,h' 

+ f dXidx2dX3 \i^hiXp)h'{Xi,X2,X3,l - Xi- X2- X3)\^ ^ . (22) 

° h,X,p,h' 

for any Pi,P2. This also ensures that the light-cone momentum sum rule is satisfied, even 
at finite DLCQ cutoff K, since translation invariance in the x~ direction is preserved by 
DLCQ. 

Since there is 90° rotational symmetry about x^ for a state with P = 0, it is possible 
to distinguish the angular momentum projections jTs mod 4. There is also exact symmetry 
under ^-parity, charge conjugation C, and transverse reflections in the x^ and x"^ direction, 
Vi, V2- Although the parity V = 'Pi'P2'P3 is dynamical and in general broken, one can 
associate a parity to boundstates from their behaviour under the free particle limit of V3. 
Indeed, there is a Z2 kinematic symmetry 

Vfiphh'il -x,x) ^ iphh' {x, l-x) , (23) 
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which corresponds to the free V3 operation in the zero-hnk sector, that is exact at any cut-off 
K. In this way, one has enough information to identify the J'^'^ structure of hght states 
unambiguously. 

In general, we will find that the 0~"^ pion mass is split from the 1 rho mass. Due 
to violations of covariance, the J^s = component of the rho (po) "will also split from its 

= ±1 components (p±) which are always degenerate on the transverse lattice at P = 0. 
In view of the low-energy nature of the truncation of the colour-dielectric expansion, we do 
not analyse heavier mesons, although their eigenfunctions are obtained as a by product of 
our calculations. 



C. Renormalisation 

We have constructed a gauge theory with transverse lattice and Tamm-Dancoff cut-offs 
that we do not intend to extrapolate and a DLCQ cut-off that we do. The first step in 
the renormalisation process is to ensure finiteness of physical observables in the limit K — >■ 
00. It turns out that divergences exist but they require only infinite and finite self-energy 
counterterms that renormalise existing parton mass terms in the light-cone Hamiltonian. 
The remaining cut-offs that are not extrapolated obviously violate Lorentz covariance. This 
violation can however be minimized by appropriate finite renormalisation all of the couplings 



appearing in (p!OD. In this section we describe our procedure for performing these finite 
and infinite renormalisations. 

It is convenient to use one of the parameters of the Hamiltonian to set the dimensionful 
scale of the theory and define dimensionless versions of the others. Conventionally we will 
use G to set the scale, which has the dimensions of mass. It will later be related to the QCD 
mass scale by calculation of the heavy source potential p. The following dimensionless 
parameters are then introduced: 

i^-U^ = 1,2,4) ; (24) 

Since we will need to study the meson eigenfunctions of P~ as a function of and P, 
let us write, for these eigenfunctions, 

2P+P- = M'^ + i?(P) . (25) 

such that -R(O) = 0. Ai^ is the invariant mass (squared). We begin with P = 0, in which 
case the non-zero Fock space matrix elements of the dimensionless invariant mass operator 

iiVu h), (1/2, 0-), ■ ■ ■ , iVn-u t), iVn, h2)\2P^P' /G'^Kxi, h[), (x2. A), ■ ■ ■ , (x„_i, p), (x„, h'^)) 

(26) 

are enumerated in figs. 0(i)-(xiii) and table |I[ A number of comments are necessary to 
explain these amplitudes. We have defined 

Rot[A,p] = e|A||p|Sgn[A]Sgn[p] . (27) 
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In the planar diagram vertices of Fig. ^ light-cone momentum fraction {x,y,z), quark 
helicity {h,h'), and link-field orientation (A,p, cr, r) labels are given where necessary. Lines 
with a bar denote the ^'''-instantaneous propagators dz^ and dz^ for v quarks and gauge 
fields respectively. 'P' denotes that a principal value prescription is used when integrating 
light-cone continuum momentum fraction across singularities. For simplicity, we have not 
shown vertices involving only anti-quarks, which are similar to those involving only quarks. 
To these diagrams we add planar spectator lines which go to make up the full gauge singlet 
Fock state. 

At finite transverse lattice spacing a, but before the light-cone DLCQ cutoff K is imposed, 
the theory is behaving like a continuum 1 + 1-dimensional gauge theory coupled to a set 



of fundamental fermion and adjoint scalar fields [|T2|. Although super- renormalisable in the 
1 + 1-dimensional sense, the light-cone quantisation in light-cone gauge introduces its own 
characterstic divergences due to the presence of non-local instantaneous interations. Those 
originating from the instantaneous gluon propagator l/d"^ are dealt with by the principal 
value prescription in the manner established by 't Hooft ||T3[. Those originating from the 



instantaneous quark propagator have been studied by Burkardt [|T4|, whose analysis 
we briefly recall. 

A basic one-loop logarithmic divergence occurs in the quark self-energy as represented 
in the Light-cone Perturbation Theory diagram of Fig. ^(i) as the quark loop momentum 
vanishes. The cubic vertices are of the same type, with coupling either rrifka or mjks, once 
the orientations of the intermediate link fields have been summed over. The divergences are 
cancelled, in these diagrams and any others obtained by adding spectators, by an infinite 
quark 'kinetic' mass counterterm in the Hamiltonian (Fig. ^(ii)) 

(kl + k^) f-dy ^28) 



TT -Jo y 

This is not sufficient for the divergences in the two- loop diagrams of Fig. §(i)(ii)(iii) to 
cancel. One may add a finite kinetic mass counterterm Sm^, adjusted at order to 
produce finite results when Fig. §(iv) is included. Higher-loop generalisations of the same 
diagrams are also rendered finite by adjusting at higher orders in ka and kg. Dressing 
loop diagrams with instantaneous gluon lines (e. g. Fig. ^ renders them individually finite. 
As in Ref . |]l4l , our own checks of these statements for the transverse lattice theory have been 
done only in perturbation theory, but we will assume they are true to all orders. It was also 
shown in Ref. fl^, by means of simple cases, that choosing the correct counterterm (5m^ was 



equivalent to restoring parity invariance, which is not manifest in light-cone co-ordinates. 

By adjusting the finite counterterm 5m^, one ensures that the K —>■ oo limit can be taken 
when DLCQ is used. However, it was pointed out by Burkardt that, while this ensures a 
finite answer for the K oo limit of the self-energy, the use of a momentum-independent 
mass counterterm in DLCQ will not yield the same as the covariant answer for the same 
couplings. In effect DLCQ produces a finite violation of covariance. This is one of a number 
of sources covariance violation in our calculation. Rather than analysing how one might 
minimize the individual violations — it is not obvious which are the most significant — 
we will perform overall covariance tests on the boundstate wavefunctions that are the end- 
product. 

A Tamm-Dancoff cut-off on the maximum number of link fields in a state also violates 
covariance. In principle, this can be compensated by introducing spectator-dependent coun- 
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terterms [0]. In practice, that will lead to too many couplings for viable calculation at a 



physically reasonable choice of Tamm-Dancoff cut-off. However, it is necessary for finiteness 
of the quark self-energy to use spectator-dependent finite counterterms 5m?. Therefore, we 
must introduce separate counterterms 5rn^ for the Fock sector containing p links (note that 
the sector with p maximum has no finite or infinite quark self-energy counterterms). These 
are adjusted to produce finite quark self energy in addition to optimization of covariance 
of hadron wavefunctions. Since we work at the level of hadrons, the quark self-energy is 
tested indirectly. A tachyonic quark self-energy, whether divergent or not, would be sig- 
nalled by tachyonic behaviour in the lightest hadron mass. Therefore, we test for absence 
of such a divergence in the pion mass as K — > oo. A positive divergent quark self-energy 
would artificially suppress the lowest Fock sectors, that are subject to loop corrections and 
counterterms, in the hadron wavefunction as -ft" — oo (the hadron mass may remain finite). 
Therefore, we test for absence of this suppression, in particular by fitting which is a 
measure of the p = Fock component. 

In summary, taking the K oo limit with a Tamm-Dancoff cut-off on link fields in place, 
one must introduce infinite and spectator-dependent finite self-energy counterterms. Even 
though the theory is now finite, Poincare covariance is still violated by the finite transverse 
lattice spacing a, the Tamm-Dancoff cutoff, and by the use of momentum-independent finite- 
mass counterterms 5rn^. We propose to minimize these violations by finitely renormalising 
all the couplings available in . 

For perfectly relativistic dispersion, -R(P) = |Pp for every eigenfunction in eq. (p5D; this 
will receive corrections on the coarse transverse lattice. To quantify the covariance violation 
we will expand the dispersion relation for each boundstate 

i?(P) =c2|P|2 + 0(P^) . (29) 

The transverse speed of light c will in general differ from one (the speed in the direction). 
A simple criterion, which worked well in previous studies, is to minimize this difference in 
the low lying eigenstates of P~, ignoring the anharmonic terms in R. 

The same procedure may be carried out for glueball boundstates to constrain the pure- 
glue interactions in the Hamiltonian, independently of the meson sector (at large N^). In 
addition, the rotational invariance of the potential between heavy sources may be optimized. 
These latter tests have been described in detail in previous work P]. We will use the string 
tension ^/a from the potential to set the QCD scale from experiment. 

The transverse lattice lagrangian (^) contains terms that also violate chiral symmetry 
explicitly, via the couplings irtf and ks- Since we work at the level of hadrons, a measure of 
chiral current non-conservation is provided by the pion mass in a covariant stable theory, as 
a result of the PCAC theorem. This measure loses accuracy if the theory also has significant 
explicit covariance violation, as in our case. The explicit violation of chiral symmetry could 
be minimized by tuning further chiral-symmetry violating counterterms, which we discuss 
in the appendix. However, these lie at higher order of the colour-dielectric expansion. In 
the present calculation, we finitely renormalize the hamiltonian to fit the experimental pion 
mass, since this is naively a measure of chiral current conservation. Since the explicit 
violation of chiral symmetry is actually larger than that suggested by M.t^, we find that we 
must also fit the experimental rho mass in order to maintain a realistic pi-rho splitting. 

Thus, in addition to optimizing covariance via boundstate dispersion, we are proposing 
to fit four experimental numbers M.-,,, M.p, fn, and y/a in order to accurately determine 
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the couplings in our effective hamiltonian P~ . Since QCD with degenerate ffavours contains 
only two fundamental parameters, the transverse lattice hamiltonian is not determined from 
first principles. However, as described above, direct tests of parity and chiral symmetry 
might allow one to reduce the number of phenomenological parameters further. We leave 
this for future work. 



III. DETERMINATION OF HAMILTONIAN PARAMETERS 

In order to reduce the number of coupling variables in the minimization process, this 
is done in two stages. First, we examine glueball eigenfunctions of P~ , that contain only 
link-fields, and the rotational invariance of the groundstate potential between two heavy 
sources of colour. Here we follow, with one exception, exactly the same procedure used 
in refs. and so omit all details. The exception is that, instead of using anti-periodic 
boundary conditions for link fields in x~ , we re-did the calculations with periodic boundary 
conditions in order to be consistent with the conditions used in the meson sector later. 
Note that these 'pure-glue' calculations extrapolate both K and the Tamm-Dancoff cut-off, 
constraining very precisely the couplings in P~ relevant to that sector (/i, I2, 13, b, mi,) when 
covariance is optimized. It is not necessary at this stage to use any phenomenological input. 

We searched for a trajectory in coupling space that optimized the Poincare covariance 
of glueball wavefunctions and the potential between heavy sources of colour. A fairly well- 
defined one-parameter trajectory is picked out. We chose the best point on that trajectory, 
which in effect fixes a, corresponding to a value of the link field mass mb = 0.276. At this 
point, we find G ~ 2.75A/a" and aG ~ 4, where a is the string tension of the asymptotically 
linear potential found between two heavy sources. If one takes ^/a = 440 MeV, then 
G ~ 1200 MeV and a ^ 2/3 fm. The values of the other couplings determined by this point 
are shown in table |I[ 

Having fixed a subset of the couplings, we fix the remaining ones sensitive only to the 
meson sector. We investigated the Tamm Dancoff-cutoff up to four links, but show results 
for a three-link cut-off, since a better sampling of couplings is achievable in this case. The 
transverse speed of light c is optimized in the dispersion of the pion and each component 
of the rho, together with the difference between the calculated mass and the physical 
value. As described in the previous section, we find we must include fits to the physical 
values of A^p and f-„ in the optimization procedure in order to accurately pin down the 
remaining undetermined couplings of the Hamiltonian, which are shown in table y. 

Table |T| shows information on the pi and rho states at these couplings. One notes that the 
spin ±1 projections of the rho still badly violate covariance, splitting the Lorentz multiplet 
and having asymmetrical dispersion. Since the rho is not yet behavingly covariantly overall, 
we do not attempt a detailed phenomenological analysis of the resulting wavefunctions. On 
the other hand, we are able to achieve a relatively symmetrical dispersion for the pion, with 
intercept A^^ = 171 MeV and decay constant = 132 MeV, close to the experimental 
values. (The exact values would be obtainable with a sufficiently fine sweep of couplings). 

We checked that no Fock sectors are being artificially suppressed and that the truncation 



to no more than three links is not causing severe 'finite-volume' effects. Table |I^ shows that 
the peak in the transverse spatial distribution of the pion wavefunction is well-accomodated 
by the three-link cut-off (results are similar for the rho). However, the tail of the wavefunc- 
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tion at four and higher hnks may contain a significant total probabihty, which will affect 
observables sensitive to very small transverse momenta. Therefore, in this paper we re- 
strict our attention to observables integrated over all available transverse momenta. Indeed, 
when varying the Tamm-Dancoff cut-off above three links, we find very little change in the 
observables investigated below. 



IV. PION OBSERVABLES 

A. Valence quark structure function 

The valence quark distribution function is defined as 

Vix) = J2\'^hh'ix,l- x)]"^ + dy \'iphix)h'ix,y,l- X -y)]"^ 

h,h' A h,h''"^ 

rl—x rl—x—y 

+ dz \^ph(xp)h'ix,y,z,l- X -y - z)]"^ ^ . (30) 

\,ph,h'-'^ -'^ 

It is the probability for a quark to carry light-cone momentum fraction x. The result we 
find for V{x) in the pion on the transverse lattice in the three-link truncation is shown in 
Fig. |[ The raw (discrete) DLCQ data for K = 10, 12, 15, 20 are displayed together with an 
extrapolation to ^ oo. To produce this, at each K data is fit to the form 

xV{x) = {l-xfx^ia + by^ + cx) . (31) 

xV is directly the momentum distribution and is easier to extrapolate since it vanishes at 
a; = We note that the simple form x°'{l — x)^, used to parameterize early experimental 
data, is not sufficient to fit our result. It is necessary to drop the x = 1/2K and x = 1 — 1/2K 
points from this analysis since they do not join smoothly to the rest of the distribution. This 
is because endpoint data are prone to artifacts resulting from the vanishing of some of the 
interactions in Table |I|. The smooth curves at each K are then extrapolated pointwise, by a 
(good) fit to a quadratic in for a large set of values in the interval 0.1 < x < 0.9. The 
grey region represents the uncertainty from the extrapolation only. 



The extrapolated data fits the form ( pl| ) with 13 = 0.33(2), a = 0.3(1), a = 0.33(3), 
h = —1.1(2), c = 2.0(3). The errors are from the extrapolation only. Bearing in mind that 
the extrapolation is based on fits to data that do not cover the endpoint regions, the true 
errors on a and (3 are likely to be much larger. From the first moment = Jq xVTr{x)dx, 

we find that 32% of meson light-cone momentum is carried by the quarks, with the same 
carried by the anti-quarks. In the range 0.1 < x < 0.9, over which there is some measure 
of control, the result for is reminiscent of the constant VJr = 1 distribution resulting 
from the chiral limit of chiral quark models (see Ref. |]TB[ and references therein). However, 
because of the rapid rise at small x, which is expected on general grounds from Regge-type 
behaviour fl^, the fiat part of the distribution is at ~ 0.7. Moreover, in the chiral quark 
models, all the light-cone momentum is carried by quarks, while here 36% is carried by the 
link fields representing gluonic degrees of freedom. 

A well-defined transverse scale is associated with V{x) above, namely, the transverse 
lattice spacing a. If we were to repeat the calculation at a different a, one would expect 
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to see an evolution of as a result of the changing wavef unctions. In practice, the current 
transverse lattice method is only able to explore a small window in a — small enough to 
suppress discretization errors but large enough for the use of massive disordered link fields 
— which is too small to reliably see evolution. Perturbative evolution equations typically 
use a different renormalisation scheme, so there is no simple matching between scales. This 
problem is common to most non-perturbative approaches to QCD, except Euclidean lattice 
QCD where (moments of) V can be calculated in schemes matched to perturbation theory 
p!8| . In low-energy effective theories for pion structure, such as QCD sum rules [jT9|, pO| , 
chiral quark, colour-dielectric transverse lattice, truncated Dyson-Schwinger etc., one 
is reduced to fixing the scale for input to perturbative evolution equations heuristically by 
demanding, say, agreement of the first moment of V with experiment at some scale. 

If we demand that (xV^) ~ 0.21 at a scale of 2 GeV, as suggested by the analysis of 
E615 and NAlO pion-nucleon Drell-Yan data by Sutton et. al. [^] for the valence quark 
distribution, then the scale associated to our result, if it were used as input for leading order 
non-singlet evolution, is /i ~ SOOMeV; this is reasonable given that a^^ = SOOMeV. In Fig. |^ 
we show our result for xV^(x) evolved to 6.6GeV and compared with the raw data for the 
valence distribution deduced by E615 by combining data over scales 4 — 8.5 GeV. For 
completeness, we also show fits to xV = a;'^(l — x)^ produced by earlier experiments NAlO 
and NA3 pSf. In the valence region x > 0.5, our result agrees with the most recent 



experiment, which claims a more accurate representation at large x. At smaller x there is not 
much agreement, either between experiments or with our result. This is hardly surprising 
given the sensitivity of this region to assumptions about the sea quarks or their measurement. 
In fact, our calculation contains no sea quarks since it is at large Nc- The recently discovered 
enhancement of initial-state interactions |2^ is also expected to be most significant at small 
X, throwing into doubt the simple connection between light-cone probabilities and the Drell- 
Yan cross-section It is obviously desirable to have data on from sources other than 
the Drell-Yan process. This is also important from the theoretical perspective, given that 
the current Dyson-Schwinger approach predicts a completely different shape for the pion 



distribution function p8 



B. Distribution amplitude 

The distribution amplitude (in A_ = gauge) for the pion is defined by 

(0| ^(z) 7^^75^1/(0) \MPn) = UP'' /'e'"(^-^V.(a:)c/a; (32) 



with the normahsation condition 



0^(x) = l. (33) 



If the quark field correlator is to be evaluated at equal light-cone time, = 0, then z = 
and z~ is arbitrary. This then measures the amplitude for zero transverse separation of 
quarks in the meson light-cone wavefunction. For the transverse lattice one finds 

?/^+_(x, 1-x) = yJ— 0^(X) (34) 
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from the 7"*" component of (p2D. Fig. ^ shows our results for the distribution amphtude at 
various K and extrapolated to iiT = cx) in a three-link truncation. The raw (discrete) DLCQ 
data has been fit at each K to the first few terms of the conformal expansion p9|, |30| 



0^(x) = 6x(l -x)[l + a^cT^^ - 2a;) + a^cT{^ - 2a;)} . (35) 

An extrapolation of the coefficients with A + E/i^ + C/i^^ yields a2 = 0.15(2), 04 = 0.04(1), 
confirming that truncation of the conformal expansion is justified. The same values are 
obtained if the fit curves at each K are pointwise extrapolated and then refit to (pSj). 
Finally, the result is insensitive to whether the x = 1/2K and x = 1 — 1/2K endpoints are 
included in the fit or not, so we have shown them in Fig. ^ also. 

The distribution amplitude is indirectly accessible through the pion transition form factor 
F-i^Y'-yiQ"^) measured at CLEO |^T|. A perturbative QCD analysis relates this to the inverse 
moment, up to radiative corrections A, 

= f' + A = 3(1 + a2 + a4) + A (36) 



Att ' Jo X 



An analysis of the data in Ref. |3^ extracted 02 + 04 = 0.05 ± 0.07 at scale 2.4 GeV, taking 
into account next-to-leading order corrections 0{as) in A. If we assume, following the 
structure function analysis, a scale 0.5GeV for the transverse lattice result, when evolved to 
2.4 GeV by the 1-loop evolution equations we find 02 = 0.07(1), 04 = 0.01(1) including only 
DLCQ errors. Although our result seems consistent with experiment, a couple of comments 
are necessary. The inverse moment is highly sensitive to the endpoint regions of (pj^, which 
are not well covered by the extrapolation of the DLCQ transverse lattice result. Also, the 
leading radiative corrections in A are large ~ 20%, so one might ask about higher order 
corrections. The reader is referred to refs. [TO], |5B| for a more detailed review of the various 



theoretical and experimental results relating to (f).„. 

Diffractive dissociation on a nucleus tt+A — > A+jets has been used to measure a cross- 
section related to (j).,^- A number a theoretical analyses of that relation have recently been 



performed |]35|, which differ in their conclusions about the precise relationship. Our result, 
when evolved to the higher transverse momentum scale of the experiments, is consistent with 
any one of the analyses, being close to the asymptotic form 6a;(l — x). We mention that 
our DLCQ transverse lattice result for 0^ is close to one previously obtained in a one-link 
truncation using very similar methods 0], although 04 was not fit and the normalisation 
Jtt was completely wrong in that case. The same 1-link truncation was investigated in 
Ref. by using basis functions instead of DLCQ and a similar (but not identical) criteria 
for fixing the Hamiltonian couplings. That gave a distribution amplitude a little closer to 
the asymptotic form, although a value for 02 was not extracted and no error estimate was 
given. Thus, we can say with some confidence that our result is neither the 'double-hump' 



first found by Chernyak and Zhitnitsky [19| using (local) sum rules nor the 'narrow-hump' 
one would deduce from the latest Euclidean lattice measurements of the lowest moment of 
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C. Quark helicity correlation 



Although the pion is spin 0, it nevertheless contains a complicated spin structure. One 
measure of this is the quark helicity correlations 



h \ h ■'^ 

rl — x rl — x—y 

+ dz \^ph(xp)hix,y,z,l- x-y - z)\'^ ^ , (37) 

X,p h ^0 ^0 

C-*i(x) = Y^\^_,,(x,l-x)\' + J2J2 r^dy\^p.,^,),ix,y,l-x-y)\' 



A h 

rl — x rl~x—y 

+ ^y dz \ip_^xp)hix,y,z,l-x-y- z)]"^ ^ . (38) 

X,p h ''^ ''^ 

They measure the probability for a quark to have light-cone momentum fraction x and 
helicity either parallel or anti-parallel to that of the anti-quark. Therefore their sum is 
normalised to one (when integrated over x). These functions are plotted in Fig. |^ for the 
pion. Although different from one another in the bulk of x, both C'^^'^''^{x) and C^'^*'(x) have 
exponents a and (3 consistent within errors with those of V^(x) (see Fig. We estimate 
that ^ 

" CP^-'^rfx ~ 0.45 . (39) 



Therefore, one is almost equally likely to find quark helicities aligned as anti-aligned! 



V. CONCLUSIONS 

We have extended coarse transverse lattice calculations to physically realistic cut-offs on 
the anti-quark — quark separation. A general light-cone hamiltonian in the large limit 
was expanded in powers of dynamical fields and we studied a truncation of that colour- 
dielectric expansion. This included all possible cubic terms and most of the quartic terms. 
By optimising Lorentz covariance of glueball, heavy-source and meson boundstates, the 
remaining freedom in the couplings in the Hamiltonian was reduced. By studying other 
symmetries, such as parity and chirality, it may be possible to constrain them further. In 
this paper, we performed a phenomenological calculation by fixing the remaining freedom in 
the couplings to best fit ^/a, M-t^, -Mp, and f^, (two of these are parameters of first-principles 
QCD). 

The lightest meson bound state has the quantum numbers of the pion and exhibits a 
reasonably covariant lightcone wavefunction. Comparing the predictions of this wavefunction 
with various experimentally measured observables for the pion, we find consistency in the 
regions insensitive to sea quarks. New observables, which in principle can be extracted from 
a higher twist analysis of experiments, follow from the multi-parton correlations in the light- 
cone wavefunction. As an example, we computed the anti-quark — quark helicity correlation, 
with somewhat surprising results. Because the tail of the pion wavefunction is still truncated 
in our calculation, we did not compute observables sensitive to small transverse momentum. 
Nevertheless, it would be interesting to look at the general features of the skewed parton 
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distributions for intermediate momentum transfers, since little hard information is available 
for these important observables. It should be straightforward to extend the calculations to 
strange mesons and heavy-light mesons. 

There are still some shortcomings in the calculation. The rho boundstate is not yet be- 
having covariantly. Our optimization of chiral symmetry could be considerably improved. 
Given the close connection of Lorentz and chiral symmetry on the lattice, we believe that 
these problems are related. In particular, higher-order terms in the colour-dielectric expan- 
sion can fulfill a dual role to improve both these symmetries. 
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Geneva College undergraduates E. M. Watson and J. Bratt for help with developing the 
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APPENDIX: Chiral symmetry 

The lattice Lagrangian explicitly breaks chiral symmetry 

^ ^ e~'^^^^ , (40) 

through the bare mass-term fif and Wilson term Kg- The standard test for this at the hadron 
level is PCAC 

{0\d,A'^\^P^) = UMl , (41) 

where is the axial current. Without knowing the precise form of A^, one can use A^jr 
to quantify the amount of explicit chiral symmetry breaking relative to other scales, such 
as the pure-QCD mass gap or the spontaneous chiral symmetry breaking scale given by the 
difference between M.^^ and masses of other light mesons. The result (|4l|) relies on exact 
Lorentz covariance, which is not present on the transverse lattice. In fact, in the calculation 
performed in this paper, the splitting of the p Lorentz multiplet is of comparable strength 
to the TT-p splitting. This suggests that explicit chiral symmetry breaking effects are larger 
than M.T, would suggest, perhaps of the same order as spontaneous chiral symmetry breaking 
effects. 

Explicit chiral symmetry breaking could in principle be tested more directly. There is a 
1 + 1-dimensional Noether 'vector' current 

j- = ^^(x)7-vl/(x) (42) 

X 

that is conserved under the equations of motion, i.e. 9aj° = 0, and a corresponding chiral 
current 

j5" = ^^(x)7'^75Vl>(x) (43) 

X 

for which we find 

d^3t = 2/^? E hFli^)^F,i^) . (44) 
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where Fh is defined in eq. (§). One might then use a matrix element such as 



(0|9„j°|^0 (45) 

to quantify exphcit chiral symmetry violation, minimizing it by finite renormalisations of 
couplings, since the vanishing of ( ^Sj) is a necessary condition for conservation of the four- 
dimensional axial current. There are a few difficulties that must be overcome before this 
would be practical however. The expression (|4^) has a normal-ordering ambiguity similar to 



the hamiltonian. Moreover, it is much more computationally expensive to perform symmetry 
tests with eigenf unctions rather than eigenvalues. On a coarse lattice, the chiral symmetry 
breaking couplings are also strongly constrained away from zero by Lorentz covariance; for 
example. Kg is needed to avoid fermion doubling 0. It would therefore be desirable to have 
further independent chiral symmetry breaking couplings in the theory to tune. 

Natural candidates are the transverse lattice versions of the Sheikholeslami-Wohlert terms 



|37|, ^a^'' F^y^i ^ that can be used in Euclidean lattice gauge theory to remove 0(a) contri- 



butions to chiral current non- conservation |]38[. They become 

*'(x)a''* (M^(x)M,(x + ar) - Ms(x)M^(x + as)) ^(x + ar + as) , (46) 

^(x)a+-F+_^(x) (47) 

While in the dimensional counting classification of Euclidean lattice operators they occur 
along with Wilson terms at dimension five, on the coarse transverse lattice their significance 
is not so obvious since they contribute to higher orders of the colour-dielectric expansion. 
Terms of the form (^Bp(^) in the Lagrangian give rise to coupled constraint equations of 
motion for non-dynamical fields. If solved order by order in dynamical fields, they give rise 
to new interactions in the gauge-fixed light-cone hamiltonian starting at order u^, v?M'^. Of 
particular interest are interactions generated at order v?M^ and u^M that fiip the helicity 
h of quarks. Such interactions carry the spontaneous chiral symmetry breaking effects in 
effective light-cone hamiltonians PS|; the nifka and kska terms performed this function in 
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(^) 


[x+y x) 1"^/^^ + ^/'^aSgn[A]((5|A|2 l^A|l) j 


iii) 




(iii) 




{iv) 


Mx+y)^ i^l^hh' + kl^hh' {hp S-Xp i/iRot[A, p]) 
+kaksS-hh' (Sgn[A](iM|A|i - d\x\2) + Sgn[p](i/i(5|p|i - S\p\2)) | 


{v) 


47r(x+2/)VP ^ ^^-^-^ + l/lRot[A, p]) 

+kaks5^hh' (Sgn[A](iM|A|i - 5\x\2) - Sgn[p](i/i(5|p|i - ^|p|2)j | 


(vi) 


'"6 

a; 


(vii) 


-Ip 2/+2 \ 


iyiii) 


\ V^a;z2/(a;+{/-2)(2-2/)2 y 


{ix) 


+^4((^(tA + <J-o-p<5-Ar|R-Ot[cr, A]|) - 6(5Ap(5o-r |Rot[cr, p]|} 


(x) 


^.^xyzix^z^y) {^^^'^-A^Ar'^Ap + /2<5.r<5-Ap5|A||.| 

+ 'l('^"-A-T<'V^"' P\\ + f^-a//Ur R"t [fT. A] ) - M_A//^<Tr R"t [fT. A] } 


(xi) 


-1 i:y-z)(2x^y^z) ^ 
»^ ^/xzy{x+y+z){z+yy 


(xii) 


-1 {x-z){2y-x-z) ^ ^ 

^ xzy{x—y+z){z+x)'^ '^^ ''" 


(xiii) 



TABLE I: Matrix elements of the dimensionless invariant mass operator 2P+P /G^ in Fock space. 
Momentum conserving delta functions are omitted for clarity. 



rub 


b 


h 


h 


k 


kg 


ka 


ruf 


5ml 


5ml 


5m2 


0.276 


0.768 


-0.169 


-0.186 


0.024 


0.420 


0.652 


0.236 


1 


-0.127 


-0.035 



TABLE II: Optimum coupling constants at a = 2/3 fm for a three-link truncation. Note that 5ml 
was swept more coarsely than the other couplings. 
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State 


Mass(MeV) 


c 


TT 


171 


1.02 


Po 


828 


0.99 


P+ 


457 


1.04 


P- 


457 


0.76 



TABLE III: Meson dispersion at the optimum couplings. 



#Links 





1 


2 


3 


Probability 


0.097(14) 


0.661(10) 


0.150(8) 


0.087(2) 



TABLE IV: Probability for finding a certain number of links in the pion. The extrapolation errors 
in brackets are from sl 1/K extrapolation. 



h'^^" \q,)W ^^1' ^ ' ^3 ^ Ap)h' (^i, ^,x,,x^) 

X, h Xi h h 

X2 A- 

+ x^X + + 

X3 p 

^^h' X3h' x,h' 

FIG. 1: Planar diagram representation of the Fock space structure of a meson boundstate. Solid 
lines represent quarks/anti-quarks, chain lines link fields. 
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FIG. 2: Planar diagram representation of the elementary amplitudes contributing to (|2^). Vertical 
barred lines are rc "'"-instantaneous interactions. 
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(i) (ii) 
FIG. 3: (i) one- loop logarithmically divergent quark self-energy (ii) logarithmically divergent mass 

insertion counterterm represented by open circle . 
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(iii) 



(iv) 

FIG. 4: (i)(iii) two- loop logarithmically divergent quark self-energies (ii) one-loop diagram with 
infinite mass insertion (iv) one-loop diagram with finite mass insertion Sm? represented by open 
box 




FIG. 5: Finite diagram with instantaneous interaction dressing 
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FIG. 6: Pion distribution function for DLCQ cutoffs K = 10, 12, 15, 20, darker data points 
meaning larger K. K ^ oo extrapolated curve lies in the shaded region. 
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FIG. 7: Pion (valence) distribution function (times x) compared to pion-nucleon Drell-Yan data. 
Solid line: transverse lattice result evolved to 6.6 GeV. Data points: E615 experiment. Short- 
dashed line: NAIO experiment fit to x°(l — x)^ form. Long-dashed line: NA3 experiment fit to 
x°(l — x)^ form. 
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FIG. 8: Pion distribution amplitude for DLCQ cutoffs K = 10, 12, 15, 20, darker data points 
meaning larger K . K ^ go extrapolated curve lies in the shaded region. 
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